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We study the structure of the eigenstates and the dynamics of a system that undergoes an excited state quantum 
phase transition (ESQPT). The analysis is performed for two-level pairing models characterized hy a [/(n + 1) 
algebraic structure. They exhibit a second order phase transition between two limiting dynamical symmetries 
represented by the U{n) and SO{n + 1) subalgebras. They are, or can be mapped onto, models of interacting 
bosons. We show that the eigenstates with energies very close to the ESQPT critical point, Sesqpt, are highly 
localized in the [/(n)-basis. Consequently, the dynamics of a system initially prepared in a i7(n)-basis vector 
with energy £ ~ T^esqpt may be extremely slow. Signatures of an ESQPT can therefore be found in the 
structures of the eigenstates and in the speed of the system evolution after a sudden quench. Our findings can 
be tested experimentally with trapped ions. 
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Introduction.- Quantum phase transitions (QPTs) occur at 
zero temperature. They correspond to an abrupt change in the 
character of the ground state of a system when a control pa¬ 
rameter passes a critical point iQl]. The subject, which perme¬ 
ates condensed matter and nuclear physics, has become one 
of the highlights of experiments with cold gases, where tran¬ 
sitions from a superfluid to a Mott insulator ^ and from a 
normal to a superradiant phase have been observed. The 
investigations are not restricted to the properties of the ground 
state, but extend also to the dynamics of systems undergoing 
QPTs. In this context, one finds studies about the quantum 
analogue of the Kibble-Zurek mechanism |01, as well as the 
relaxation time il-li], revivals ||1,[I3], and temporal fluctua¬ 
tions at critical points. 

Recently, the concept of ground state QPT has been gen¬ 
eralized to encompass also QPTs occurring at excited states. 
These so-called ESQPT refer to a singularity in the energy 
spectrum caused by the clustering of excited levels at a critical 
energy cini]. This critical point can be reached either for a 
constant excited energy by varying the control parameter(s) 
or by fixing the latter and increasing the energy. ESQPTs 
have been investigated for a broad class of many-body quan¬ 
tum systems, such as the Lipkin-Meshkov-Glick (LMG) ifp^ 
, the molecular vibron th e nuclear intera cting bo - 

son m, the Javnes-Cummings II 21 L 1^ . the Dicke lEFlll, 
and the kicked-top models. Experimentall y, s ignatures 
of ESQPTs were found in molecular systems 11251 - 1^ . su¬ 
perconducting microwave billiards 1^ . and spinor conden¬ 
sates i^ . 

In terms of dynamics, it has been shown that an ESQPT 
leads to random oscillations of the survival probability in iso¬ 
lated systems iH, to singularities in the evolution of ob¬ 
servables and to maximal decoherence in open sys¬ 

tems lUM 1^ . Despite these works, studies of the effects of 
ESQPTs on systems’ evolutions are still scarce. 

In this Rapid Communication, we provide insights into the 


dynamics of an isolated many-body quantum system that un¬ 
dergoes an ESQPT. The system is prepared in an eigenstate of 
an initial Hamiltonian fLi. The evolution starts after the sud¬ 
den quench of a control parameter f that changes fii into a 
new final Hamiltonian described by a C/(n + 1) algebraic 
structure, 

= (1 -+ ^'HsO(n+l)- (1) 

This model exhibits a second order ground state QPT at 
= 0.2, which occurs between the dynamical symmetries 
(DSs) represented by the U{n) and the SO{n + 1) subalge¬ 
bras 10. It also displays an ESQPT at an energy i^ESOPTCO 
for Cesqpt The model represents systems of inter¬ 

acting bosons. It is built upon two types of bosons, a scalar 
and a non-scalar one. 'Hu{n) is the number operator of the 

non-scalar boson, while TLsoin+i) is a two-body (pairing) 
operator built from the second order invariant operator of the 
SO{n -f 1) subalgebra, which is rescaled by the system size, 
N. Such models have been successfully applied to problems 
of hadronic 1^ . nuclear ll^ . and molecular physics 0. 
Eor n = 1, they coincide with the LMG model 11^ in the 
bosonic form. 

We assume that the initial Hamiltonian corresponds to one 
of the two limits ofT-Lp, with ^ = 0 or ^ = 1. Thus, the initial 
state ItkjO)) is either a basis-vector |(/)( 7 (n)) associated with 
the DS U{n+ 1) D U (n) D ..., which defines the so-called 
“spherical” or “symmetrical” phase, or it is a basis-vector 
\4>so(n+i)) associated with the DS U{n + \) D SO{n + l) D 
..., which corresponds to the “deformed” or “broken symme¬ 
try” phase 0. We study the evolution of initial states with 
different values of the energy 

£ = (^'(0)|?iF|«'(0)), (2) 

and show that the rate at which these states change in time can 
be anticipated from the eigenstates structures. 



2 


The initial state |4'(0)) = \4'u(n)) with energy £ closest to 
^'ESQPT corresponds to the ground state of the U (n) Hamilto¬ 
nian. The evolution of this state is extremely slow. This hap¬ 
pens because the main contributions to its dynamics stem from 
very few eigenstates of TLp whose energies are exceedingly 
close to the separatrix that marks the ESQPT. These eigen¬ 
states are strongly localized in the ground state of the U(n) 
Hamiltonian 

Differently from ground state QPTs, the characterization of 
the different phases in ESQPTs is hindered by the fact that the 
order parameter does not vanish above or below the critical 
point. However, degeneracy patterns with respect to angu¬ 
lar momentum llfM 1^ and the structures of the eigenstates 
reveal information about the phases ifisll . Eor instance, the 
distribution of the 17(n)-components of the eigenstates below 
and above the separatrix resemble, respectively, the distribu¬ 
tions of the U (n) components of the eigenstates of 
(eigenstates of the deformed phase [SO(n + 1) DS] ) and of 
the eigenstates of TLp ^ 'Hu{n) for ? 0 (eigenstates close 

to the spherical phase [U{n) DS]) (see details in Ref. id). 
The aforementioned excited eigenstates at the separatrix are 
the ones that mark the change of character from one symme¬ 
try to the other, being highly localized in the ground state of 
the spherical configuration. In contrast, when written in the 
SO{n + l)-basis, the eigenstates close to the separatrix are 
delocalized and do not show particular structures. As a result, 
the evolution of IT'(O)) = |<?5<so(n-i-i)) is not much affected by 
the ESQPT. 

The U (n)-basis plays a special role in the study of ESQPTs. 
As we show below, the point of the transition is clearly re¬ 
vealed from the analysis of the structure of the eigenstates in 
this basis as well as from the dynamics starting from a U(n)- 
basis vector. 

Model - The findings described in this work were numer¬ 
ically confirmed for one-, two-, and three-dimensional vi- 
bron models characterized respectively by the 

(7(2), (7(3), and U (4) algebraic structures, as well as for the 
(7(2) LMG model 11^ . The results were equivalent, so we 
chose the 3D case to present the illustrations below. 

The (7(4) vibron model provides an algebraic framework 
for the full rovibrational spectrum of diatomic molecules in. 
It has two rotationally invariant DSs: (7(4) D (7(3) D SO(3) 
and (7(4) D S'(9(4) D SO{3). The (7(3) limit describes 
the vibrational spectrum of non-rigid molecules, while the 
5'0(4) limit corresponds to rigid molecules with vanishing 
vibrational transitions. The ground state QPT of this model 
was studied in ll43l] . 

The (7(4) dynamical algebra generators are bilinear prod¬ 
ucts of creation and annihilation operators of the scalar s 
and the vector boson operators, with p = 0, ±1. The 
scaled Hamiltonian, 'HpiA) = (1 ~ built 

from the number operator, n = v/3[p^ X which is 

the first order Casimir operator of the (7(3) subalgebra, and 
the pairing operator, P = N{N -b 1) — W^, that contains 
the second order Casimir operator of the S'0(4) subalge¬ 
bra; + L^, where Lfj, = '/2[p'^ x and 

Dfi = i[p^ X s -b s'( X IS- The (7(3)-basis, 


\[N]nL), has quantum numbers n = 0,1,..., — 1, 

and L = 0 or 1,..., n — 2, n. The 5'(7(4)-basis, | [N]'wL), 
has quantum numbers w = 0 or 1,..., N — 2, N and L = 
0,1,... ,w — l,w. The matrix elements of the two Casimir 
operators in any of the two bases can be found in 
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FIG. 1: (Color online) Top left: normalized excitation energies V5 G 
N = 200. Top right: density of states p; N = 400, ^ = 0.6. Middle 
panels: participation ratio of the eigenstates in the U (3)-basis; N = 
600 = 0.6 has also N — 2000 (bottom curve) and a zoomed in 

inset]. Bottom panels: same for the 5'0(4)-basis. All panels: L — 0; 
vertical lines mark Eesqpt from Eq. (1^. 

Separatrix and density of states.- Eor ^ above the criti¬ 
cal point, = 0.2, there appears an energy region in the 
spectrum of 77(7(4) with a large density of excited levels 
signalizing the ESQPT. This is illustrated in the top left panel 
of Eig. □ where we plot the normalized excitation energies 
Ek /N for all levels as a function of the control parameter Q 
[Throughout this work, the value of the energy of the ground 
state is set to zero. In the figures, the units for energy and 
time are arbitrary.] The dashed line is the separatrix that di¬ 
vides the states with different physical characters; those closer 
to the deformed configuration below the separatrix and those 
closer to the spherical configuration above it. The equation 
for the separatrix as a function of ^ in the mean field (large 
N) limit is 

T^esopt = (1 — 5^esqpt)^/(16^esqpt)- (3) 
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This equation was derived in for the C/(3) model and 

can be extended to other U(n) models ll^ . The large density 
of states along the separatrix is made evident by the peak at 
i^ESQPT in the energy levels histogram in the top right panel of 

Fig.IU 

Structure of the eigenstates - To analyze of the structure of 
the eigenstates of the total Hamiltonian [/( 4 )> we use the par¬ 
ticipation ratio P, which quantifies the level of delocalization 
of a state in a particularly chosen basis lli^ . A large value 
indicates an extended state in that basis and a small value, a 
localized state. For an eigenstate written in the C/(3)-basis, 

m = EtLcL^^mnL),, 
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For eigenstates written in the S'0(4)-basis, I'fk) = 
|[A^]u’T)fc, we have Note that the sums 

in n and w are in increments of two units. 

When ^ = 0, the eigenstates coincide with the U (3)-basis 

(k) 

vectors and ^’( 7 ( 3 ) =1- As ^ increases, the average level of 
delocalization in the ?7(3)-basis grows. However, the depen- 
dence of the values of ^*^ 3 ) on the energies Ek changes sig¬ 
nificantly for ^ before and after the critical point. This can 
be seen in the five middle panels in Fig. [T] where we plot 
Vi Ek/N. Across the region 0 < ^ < ^ 

smooth function of energy. In contrast, a singularity appears 
above the critical point; a pronounced dip becomes noticeable 
at i^ESQPT- Its location moves as ^ increases, following the 
increasing value of E'esqpt- As seen for ^ = 0.6, the dip be¬ 
comes more pronounced as the value of N increases. 

The bottom panels of Fig.[T]depict P^o{i)/^ vs Ek/N for 
different values of The average level of delocalization in 
the 5'0(4)-basis decreases as ^ increases and = 1 

when ^ = 1. At ^ = 0, the plot shows a dip in the middle of 
the spectrum with no relation to the ESQPT. Below the critical 
point, as ^ increases, this dip moves toward lower energies and 
fades away. Above a discontinuity appears at i^ESQPT- It is 
much less conspicuous than the dip in the U (3)-basis and it re¬ 
quires large N to be apparent, as seen in the inset for ^ = 0.6. 
Notice also that PgQ^^'j peaks to its maximum value for the 
states with energies right above Eesqpt- From these observa¬ 
tions, one concludes that PR in the two bases can be used to 
identify the ESQPT point, but it is by far more evident in the 
C/(3)-basis, especially for small N. 

At first sight, the drop in the value of Pu{ 3 ) and Pso{a) 
at E^esqpt is counterintuitive. PR usually reflects the density 
of states. Even though the values of PR are intrinsically at¬ 
tached to a basis, in general, one expects the states to be more 
extended where the density of states is larger. To better un¬ 
derstand what happens to the structure of the eigenstates in 
the vicinity of E’esqpt, we analyze in Pig. |2] for ^ = 0.6 and 
N = 600, the contributions to |'0fc) from each basis vectors 
for three eigenstates, one before the ESQPT {k = 49), the 
eigenstate closest to the ESQPT (fc = 148), and one after the 
ESQPT (fc = 249). We depict the squared coefficients p 


(middle panels) and \Cw'^ p (bottom panels) as a function of 
the energies of the bases, ej, = {[N\hL\%-njjMN]hL) with 
b = n, w. A similar study was done in Ref. If^ . but the plots 
were for the coefficients with respect to the index of the basis 
vectors, instead of their energies. 



Normalized SO(4) Basis State Energy (e^yN) 


FIG. 2: (Color online) Top panel: ratio R vs N for the eigenstate 
in the (7(3)-basis closest to the separatrix (top curve) and for the 
one with the largest Pu( 3 ) (bottom curve). Middle and bottom pan¬ 
els: squared coefficients of the eigenstates, respectively, in the U (3)- 
and SO(4)-basis vj the energies of the corresponding basis vectors; 
N = 600. All panels: L = 0 and ^ = 0.6. Left: low energy 
eigenstate k — 49, Ek/N = 0.1767. Middle: eigenstate with en¬ 
ergy closest to the separatrix, k = 148, 17esqpt/N = 0.4173. Right: 
high energy eigenstate k = 249, Ek/N = 0.6513. Vertical lines 
mark the ESQPT critical energy from Eq. 0. 

The structures of three selected eigenstates in the (7(3)- 
basis are illustrated in the middle panels of Pig. |2l Por an 
eigenstate with low energy (fc = 49), the largest contributions 
(largest come from low values of e^, with a peak 

at each edge of the contributing energy interval. As Ek in¬ 
creases and approaches Eesqpt, the largest value of 
on the left of the contributing energy interval moves away 
from the boundary, towards higher energies, and its ampli¬ 
tude increases substantially. At the ESQPT {k = 148), this 
enhanced peak obscures the other components. Compare the 
y-axis scale for k = 148 with k = 49, 249. Pinally, as Ek 
further increases beyond the critical point, the enhanced peak 
decreases and moves back to the left boundary of the con¬ 
tributing energy interval, the whole interval naturally moving 
towards larger values of e„, as seen for k = 249. 

The eigenstates in the U (3)-basis with energy very close 
to the separatrix have a blunt preference for the first basis 
state, whose energy e„=o N also ^ Eesqpt- This is em¬ 
phasized by the top panel of Pig. |2l which shows the ratio 
R = \Cn=o/Cn=v\'^ (where \Cn=v\^ is the second largest 
component of the eigenstate) for various system sizes. R is al¬ 
ways larger than 1 for an eigenstate close to the separatrix (top 
curve). In contrast, i? ~ 1 for eigenstates with large PR (bot¬ 
tom curve). The favoritism for |[A^]0, 0) can be understood 
as follows. At the separatrix, the structure of the eigenstates 
change from being close to the 50(4) deformed symmetry to 
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being more spherical [[/(3) DS]. The eigenstate IV’esqpt) at 
the separatrix is the ground state of the spherical {U (3)] con¬ 
figuration, so its largest contribution comes from e„=o- The 
high level of localization of |'i/'esqpt) in the C/(3)-basis can 
also be explained with dynamical considerations based on the 
classical limit of the model lld^ . 

The structures of the eigenstates in the S'0(4)-basis are less 
striking (bottom panels of Fig.|2]i. For an eigenstate with low 
energy {k = 49), the largest amplitudes of Cw occur at low 
values of tw As Ek increases, the contributing energy inter¬ 
val moves towards larger values of Cu,. 

Dynamics.- We start by analyzing the effects of the ESQPT 
on the time evolution of initial states that correspond to U (3)- 
basis vectors, given the special role of this basis, as described 
above. Thus, T-Li = ’H( 7 ( 3 ) TLp = 'Hu(^ 4 ). The simplest 
quantity to evaluate how fast an initial state |T'(0)) changes in 
time is the survival probability. 


Sp{t) = 




k 


2 


(5) 

where = {ipkliNjnL). V(t) is the discrete Fourier trans¬ 
form in energy of the components Thus, the distri¬ 
bution of Ek weighted by | | ^ for a chosen initial state 

characterizes the decay of the survival probability lld^ . 

In Fig. [3 (a), we show the weighted energy distribu¬ 
tion for an initial state with n = 0. Its energy, £n=o = 
([7V]0,0|?f[F|[A^]0, 0), is the closest one to E'esqpt- As ex¬ 
pected from the analysis of the structure of IV'esqpt), the dis¬ 
tribution is strongly localized at i^ESQPT- Thus, the evolution 
of this first U (3)-basis vector is very slow, as seen in Fi g. [3 (d) 
and it does not accelerate as the system size increases foil . 



FIG. 3: (Color online) Weighted distribution of normalized excita¬ 
tion energies for three [/(3)-basis vectors: n = 0 (a), n = 2 (b), and 
n — 1332 (c). In (d): Survival probability for the three basis vectors 
evolving according to Hu( 4 ) ■ A = 2000, L = 0, and ^ = 0.6. 

The second C/(3)-basis vectors for L = 0, |[A]2,0), is 
slightly more spread out in energy [Fig. |3](b)] than | [A]0,0), 
its dynamics being then faster [Fig.|3(d)]. As n increases, the 
energy distribution stretches further in the direction of ener¬ 


gies smaller and also larger than Aesqpt- As a consequence, 
the second state with closest to Aesqpt is not |[A]2,0), 
but instead a high-n basis vector with a more or less symmet¬ 
ric distribution around Aesqpt [Fig-|3](c)]. Its decay is much 
faster than that of the first basis vectors [Fig. |3 (d)] and, for 
the system sizes studied, it gets faster as N increases B. 

Similarly to what happens to the structure of the eigenstates 
in the S'0(4)-basis, the projection of one \ [N]wL) onto {ipk) 
has an ‘accordion-like’ behavior as increases from zero. 
The energy distribution of E^ weighted by p is localized 
around low (high) values of Ek when Sw is small (large) and 
it spreads out for away from the edges of the spectrum. 
The dynamics reflects these distributions, being slower and 
the fluctuations after saturation being larger for initial states 
with energies closer to the border of the spectrum than for 
those away from it. No special behavior is observed for ^ 
Aesqpt- 

Experimental realization.- Recent experiments with 
trapped ions studied the quench dynamics of systems 

where the range of the interaction was tunable. One of 
the systems considered was described by the Ising spin-1/2 
Hamiltonian with a transverse field, which in the limit of 
infinite-range interaction corresponds to the LMG model |[l3 - 
[H. Thus, the experimental setup to compare the speed of the 
evolution for different basis vectors is already available. 

Conclusion.- In general, the dynamics of initial states with 
energies very close to the edges of the spectrum is much 
slower than for states with energies closer to the middle of 
the spectrum lld^ . Here, we showed that in an isolated sys¬ 
tem undergoing an ESQPT, a slow time evolution can occur 
also for initial states with energy very close to Aesopt- This 
is the case of an initial state corresponding to the ground state 
of the U (n) subalgebra of a model described with a (7(n -F 1) 
algebraic structure. This behavior reflects the structures of the 
eigenstates of the U{n + 1) Hamiltonian close to the sepa¬ 
ratrix, which are highly localized in the ground state of the 
U{n) Hamiltonian. Our findings have therefore identified two 
additional methods to detect the presence of an ESQPT, by 
analyzing the level of delocalization of the eigenstates in the 
U (n)-basis and by comparing the speed of the evolution of 
different U (n)-basis vectors evolving under U{n + 1) Hamil¬ 
tonians, the latter being more accessible experimentally. 

We plan to extend the present studies to consider also the 
effects of an ESQPT associated with first order phase tran¬ 
sitions SB- We also intend to study the influences of 
ESQPTs on coupled systems lIB . 
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